Intercellular synchronization of diffusively coupled astrocytes 



Md. Jahoor Alam^, Latika Bhayana^, Gurumayum Reenaroy Devi^, 
Heisnam Dinachandra Singh^, R.K. Brojen Singh^, B. Indrajit Sharma^. 

^ Centre for Interdisciplinary Research in Basic Sciences, 
Jamia Millia Islamia, New Delhi 110025, India 
^Department of Physics, Assam University, 
Silchar 788 Oil, Assam, India 

(Dated: December 16, 2010) 
ABSTRACT 

We examine the synchrony of the dynamics of locahzed [Ca^'^]i oscillations in internal pool of 
astrocytes via diffusing coupling of a network of such cells in a certain topology where cytosolic 
Ca^^ and inositol 1,4,5-triphosphate (IPS) are coupling molecules; and possible long range inter- 
action among the cells. Our numerical results claim that the cells exhibit fairly well coordinated 
behaviour through this coupling mechanism. It is also seen in the results that as the number of 
coupling molecular species is increased, the rate of synchrony is also increased correspondingly. 
Apart from the topology of the cells taken, as the number of coupled cells around any one of the 
cells in the system is increased, the cell process information faster. 

.ftTEYWORDS : Cell signaling, synchronization, diffusive coupling, network topology, chemi- 
cal coupling 
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I. INTRODUCTION 



The astrocytes in the central nervous system have var- 
ious important roles, namely, taking active part in signal 
processing [H-Q , interact with the neighbouring neurons 
[3i6| etc. which leads to important responsibility of the 
cells in predicting disease states Since these cells 
interact with the environment, Ca?^ waves propagated 
through the cells display oscillation in the nonpropagat- 
ing internal stored [Ca?'~^]i inside the cells H. These 
oscillations are sustained due to interaction of inositol 
1,4,5-triphosphate (IP3) with extracellular, cytosolic and 
endoplasmic Ca?'^ through inositol cross coup ling and 
calcium induced calcium release mechanisms |9|, llCf . 

The coupling among a network of these astrocytes in 
a certain topology is done through the process of cell 
signaling. This cell signaling is being considered as a 
means of complex communication and information pro- 
cessing among individual astrocytes, astrocytes with neu- 
rons, governing basic cellular activities and to coordinate 
various actions involving various complex coupling mech- 
It has been predicted that Ga^^ variation in 
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anisms 

IPS concentration is necessary for the intercellular Ca^+ 
waves propagation and this variation is initiated by IPS 
diffusion through gap junctions to communicate neigh- 
bouring cells |12| . Therefore synchrony in the dynamics 
of the local [Ca^+J^ oscillations is believed to be due to 
chemical coupling i.e. due to exchange of cytosol Ca?^ 



waves and IPS [iMll and can exhibit synchronization 
of the cells over long distances [l^. However, this idea 
of synchrony of [Ca?^]i oscillation by chemical coupling 
is considered to be waek by claiming that these chemi- 
cal wave propagate relatively slow and coupling through 
release of Ca?'^ is very weak as the effective diffusion of 
it is limited to very short distance [13] ■ They proposed 
electrochemical rather than chemical coupling, in which 
strong electrical coupling combined with weak chemical 
coupling, is responsible for the synchrony of these cells 
and is an effective means of long range signaling [T7| . 
However, it has been predicted that in astrocytes inter- 
cellular Ca^"*" waves can travel over several hundred mi- 
crometers and may prodive a long range synchrony [l[. 



The aim of this work is to try to raise this issue and try 
to look for a reasonable solution regarding synchrony of 
astrocytes via chemical coupling and possible long range 
interaction among the cells. We organize our work as 
following. We first study the model of Ca^+ oscillations 
in astrocytes developed by Houart et. al. [l^ in sec- 
tion II. We introduce possible diffusive coupling mech- 
anisms among a topological network of coupled cells to 
see whether there is synchrony behaviours in the local 
[Ca?'~^]i oscilations and look for long range communica- 
tion among them. We present our simulation results 
based on the model discussed in section HI and some 
conclusions based on our results are drawn in section IV. 
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II. MATERIALS AND METHODS 

The basic single cell model which is described in Fig. 
1 involves three key variables; the free calcium concen- 
tration in the cytosol (X), the concentration of stored 
Ca^^ in the internal pool (Y), and the inositol 1,4,5- 
triphosphatc, IP3 (Z) [li-l2l|. The time evolution of 
these three variables is governed by the following differ- 
ential equations: 



dX 

~dt 
dY 

~dt 
dZ 

H 

V5 
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where, Vq indicates the constant input of Ca^^ and Vi 
returns to maximum rate of stimulus-induced influx of 
Ca^'^ from extracellular medium. The parameter /3 is the 
degree of stimulation of the cell, Vi and V2 are the pump- 
ing and release of Ca^^ from cytosol to internal store and 
internal store to cytosol respectively in the Ca^^ induced 
Ca^'^ release (CICR) process with the maximum rates 
Vm2 and Vm3- K2, K^, Ky and are the threshold 
values release, pumping, activation of release by Ca^"*" 
and IPS (Z) respectively, kf \s the passive, linear leak 
rate constant of X and Y; k is the rate of Ca?^ diffuse 
into extracellular medium; V4 relates to rate of stimulus- 
induced synthesis of Z; V5 is the phosphorylation rate of 
Z by the 3-kinase; Vm2 is the maximum value and half 
saturation constant ; Kd corresponds to the threshold 
Ca^^ level; e is the reflected Z to mobilize in Ca^^; m, 
n and p are Hill co-efficients. 

The oscillations in in the variables X, Y and Z ex- 
hibit various types, namely simple oscillation, bursting, 
chaotic and quasiperiodic subject to different values of re- 
action constants and parameters in the set of equations 
n 0, HI] • These complex oscillations of the variables are 
in fact significantly stimulated by the variation of popu- 
lations of species IPS (Z) and stored 00^+ in the internal 
pool (Y) supported by some experimental reports [U [2^. 

Consider a group or network of TV such identical cells 
which are coupled by the exchange of free cytocelic 
Ca2+ ions (X) and IPS (Z) by diffusing through the 
ion-channels on each cell surface. Among these cells, 




FIG. 1: (A) The schematic diagram of reaction network of 
molecular mechanisms of the model of Ca^+ of the model. 
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FIG. 2: (A) The schematic diagram of diffusive coupling of a 
chain of cells via (A) Diffusing molecule, Ca^^., (B) Diffusing 
molecules, Ca^^ . and IPS. 



the information of individual cell reactions is transmit- 
ted via these diffusing molecules to each neighbouring 
cells in the network and a unique globally co-ordinated 
behaviour of the localized molecules i.e. [Ca2"'']i stored 
in internal pool (Y) of each individual cells in the net- 
work will be exhibited. If we consider a simple topo- 
logical network of these cells diffusively coupled in se- 
ries in the form a chain as shown in Fig. 2 (A) and 
(B); we have the following two mechanisms of diffusing 
coupling in the topology; (i) Single molecule diffu- 
sive coupling: In this case of single molecular species 
coupling mechanism, where X being taken as coupling 
molecule, the species Xi and Xi±i can interconvert via 



additional reaction channels, Xj 



X,, 



X,.. 



X,. 
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FIG. 3: The plot showing the time evolution of concentration 
of X, Y and Z showing (i) simple oscillation for the parameters 
[3: Vo=2, Vi = 2, P ^ 0.5, Vm2 = 6, fca = 0.1, Vms = 20, 
= 0.5, Ky = 0.2, = 0.2, Vms = 5, = 1, ka = 0.4, 
fc/ = 1, fc = 10, e = 0.1, 14 = 2, m = 2, p = 2 and n = 4. (ii) 
bursting for the parameters; Vb = 2, Vi = 2, /? = 0.46, Vm2 ~ 
6, k2 = 0.1, Vm3 = 20, Kxz = 0.3, Ky = 0.2, A"^ = 0.1, 
Vm5 = 30, fcg = 1, fcd = 0.6, fc/ = 1, = 10, e = 0.1, V4 = 2.5, 
m = 4, p = 1 and n = 2. (iii) chaos for the parameters: 
Vo = 2, \/i = 2, /3 = 0.65, Vm2 = 6, fca = 0.1, Vms = 30, 
/S'^ = 0.6, Ky = 0.3, ATz = 0.1, Vm5 = 50, ks = 0.3194, 
fed = 1, kf = 1, fc = 10, e == 13, 1/4 = 3, m=2, p=l and 
n = 4 and (iv) quasiperiodicity for the parameters: Vo ~ 2, 
Vi^2, p = 0.51, Vm2 = 6, = 0.1, Vms = 20, K^ = 0.5, 
Ky = 0.2, if, = 0.2, 14/5 = 30, fcs = 0.3, ka = 0.5, fc/ = 1, 
A; = 10, e == 0.1, Va = 5, m=2, p=2 and n = 4 



Xi ^i+i, Xi^i Xi, where the diffusing rates in all 
these additional reactions are k^i, kx2, kxs and kxi- De- 
terministically this corresponds to two bidirectional dif- 
fusive couplings i.e. [kxi{Xi — Xi-i) and kx2{Xi-i — Xi)] 
and [kx3{Xi — Xi+i) and kx4{Xi+i —Xi)] are respectively 
incorporated bidirectionally. So the synchronization in 
other variables Yk,{k — 1, 2, N) occurs when the rates 
kxi, kx2, kx3 and kxA are sufficiently large. Now, taking 
kxi ^ kx2 = kx3 = kxA = kx for simplicity, wc have the 
following diffusively coupled differential equations of the 
cells in the network topology we considered. 



dX, 
dt 

dt 
dZi 

~dt 



+2j-l 

j=0 

H, 



X^) 



(2) 



and (ii) Global diffusive coupling: In this case, two 
or more molecular species are considered as diffusively 



coupling molecules in the cells in the network. For ev- 
ery diffusing molecular species, four additional reaction 
channels are to be added with different rates. In the 
model we considered, we have two diffusively coupling 
molecules i.e. X and Z. Therefore we will have eight 
additional reaction channels; four for X and four for Z 
respectively. Taking the diffusive rates of each molecular 
species to be the same, the coupled differential equations 
arc given by. 



dX, 
dt 

dt 

dZ^ 
dt 



Fi + kx{Xij^2j-i — Xi) 



= G, 



H,+Y,k.{Zr+23-l~Z,) 
3=0 



(3) 



where, Xi, Yi and Zi are the variables of the ith cell in the 
chain of cells we considered. The functions in the above 
equations arc defined by, Fi = Vq + ViI3 — V2 + V3 + kfYi — 
kXi, G^ = V2-V3~ kfYi and H, = pVi - V5 - tZi. The 
parameters kx and k^ are the coupling constants of X 
and Z molecular species and arc not necessary to have 
the same value. For the topology of the finite chain of 
cells, the coupling terms containing Xq, Zq in the first cell 
and Xm+1, Zm+i in the Nth cell are neglected since these 
molecular species are not in the domain of the system we 
considered. However, if the chain become a ring by con- 
necting the two ends, then we have to apply the bound- 
ary conditions i.e. Xq — Xn, Zq = Zn and Xn-^-i = Xi, 
Zn+1 = Z\. 

The measure of synchronization of the time evolution 
of two independent and identical systems can be possible 
[2^ 1 by defining an instantaneous phase for an arbitrary 
signal r\{t) via the Hilbert transform [2^ 



J -00 t-'T 



(4) 



where P.V. denotes the Cauchy principal value. The in- 
stantaneous phase and amplitude A{t) of a given 
arbitrary signal can be obtained through the relation. 



A{t)~ 



vit) 



[{'i,j);hj = 1,2, 



- ifi{t). For any given pair of signals 
■ ,N,i ^ 7], one can therefore obtain 



the instantaneous phases (pi and 0^; phase synchroniza- 
tion is then the condition that A0 = m(f>i — ncjjj is con- 
stant with m and n being integers. Starting with differ- 
ent initial configurations, the temporal dynamics of the 
uncoupled oscillators will be uncorrelated; however upon 
coupling, the dynamics can show phase synchrony |25l - 

S3- 

Another way to measure the rate of synchrony of two 
coupled oscillators is to plot the two corresponding vari- 
ables X, x' from the two oscillators along the two axes of 
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FIG. 4: Plots of Y (left panels) and phase difference A(j> (right panels) as a function of time for different values of ke'. (i) simple 
oscillation in (a) and (b) for ke=0, 5 and 40; (ii) bursting in (c) and (d) for ke—O, 10 and 28; (iii) chaos in (e) and (f) for 
ke—O, 8 and 30; (iv) quasiperiodic in (g) and (h) for k^—O, 5 and 20; for 10 cells out of 50 cells showing desynchronized, weakly 
synchronized and strongly synchronized regimes. 



the two dimensional cartesian plane (Pecora-caroll type) 
psf . If the oscillators are uncoupled then the points in 
the plane scattered away from the diagonal. However, if 
the oscillators are coupled then the points concentrated 
towards the diagonal. The rate of concentration of the 
points towards the diagonal measures the rate of syn- 
chrony. 



III. RESULTS 

We present first one time temporal dynamics of the 
concentrations of the variables of single cell by standard 
numerical integration technique [2^ of the set of differen- 
tial equations ([Ij for various values of four sets of param- 
eters [19| given in the figure captions of Fig. 3 showing 



simple oscillation, bursting, chaos and quasiperiodicity. 
The time is measured in minutes. Then we take A'^ = 50 
identical cells, diffusive coupling is employed with X as 
diffusing coupler molecule at different rate constants at 
different times and solved the sets of coupled differential 
equations ^ as shown in Fig. 4. The results of sim- 
ple oscillations are shown in Fig. 4 (a) and (b); the 
figure (a) shows plots of the dynamics of concentrations 
of the variable Yi{i = 1,2,..., N) for the first = 10 
cells and figure (b) shows the corresponding phase plot 
(A0 verses time) for different rate constants i.e. kg — 
(uncoupled), kg = 5 and kg = 40 switched on in time in- 
tervals (0-15), (15-30) and (30-45) minutes respectively. 
The uncoupled behaviours of the cells are shown during 
time interval (0-15) minutes proved by random fluctua- 
tion of A(j) with time during this time interval and pecora- 




caroll plot of Yi verses Yio in Fig. 6 (a) where the points 
spread randomly away from the diagonal. During time 
interval (15-30) minutes the behaviours of the cells start 
co-ordinating but is still weak indicating weak synchro- 
nization which is again proved by phase plot where the 
fluctuation is reduced enormousely and by pecora-caroU 
plot whose points start concentrated enormousely along 
the diagonal. However, during time interval (30-45), the 
cells exhibit strong co-ordinated nature i.e. synchronized 
state proved by the negligible fluctuation in the phase 
plot and diagonally aligned points in the pecora-caroU 
plot. We also noticed from the results that the dynamics 
of the variables of the cells do not synchronized immedi- 
ately but takes some time after the coupling is switched 
on. 

Similarly, for the bursting case, the curves in Fig. 4 
(c) and (d); and Fig. 6 (b) present how the dynamics of 
the variable Y start co-ordinating their behaviours from 



uncoupled to weak, then to strong showing desynchro- 
nized, weak and strong synchronization. As explained 
above, the claim is based on the fluctuation rates in the 
phase plot (Fig 4. (d)) and spreading rate of the points 
from the diagonal in pecora-caroU plot (Fig. 6 (b)). The 
diffusing rates for desynchronization, weak and strong 
synchronization in this case are found to be ke = 0, 10 
and 28. In the same way, the results for chaos and 
quasiperiodic are shown in Fig. 4 (e) and (f) with Fig. 
6 (c) (for chaos) and Fig. 4 (g) and (h) with Fig. 6 (d) 
(for quasiperiodic) respectively. The diffusing rates are 
found to be different as fcg = 0, 8, 30 and fee = 0, 5, 20 for 
chaos and quasiperiodic respectively. 

Now the simulation results for the cases of simple oscil- 
lation, bursting, chaos and quasiperiodic, when two dif- 
fusively coupling molecular species namely X and Z are 
considered in the same topology of the cells, are shown 
in Fig. 5 [(a), (b)], [(c), (d)], [(c), (f)] and [(g), (h)] with 
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FIG. 6: The pecora-caroU type plot of variable Ys of dif- 
fusively coupled first and 10th cells for (a) oscillation, (b) 
bursting, (c) chaos and (d) quasiperiodic. 



Fig. 7 (a), (b), (c) and (d) respectively. The results are 
obtained by solving the set of coupled differential equa- 
tions ([3]). Similar behaviours of desynchronized, weakly 
synchronized and strongly synchronized are found as in 
the case where single molecular species diffusive coupling 
is employed. However, interestingly the mentioned be- 
haviours are found in significantly lower diffusing rates 
i.e. ke = 0,1,30 for simple oscillations, fcg = 0,5,20 
for bursting, ke = 0,5,25 for chaos and fcg = 0, 1.2, 15 
for quasiperiodic respectively. 

The results interestingly give the evidence that as 
the number of coupling molecular species are increased, 
the rate of synchronization also increases significantly. 
In other words, as the number of information carrying 
molecules (same molecular species) or (and) molecular 
species increases, the rate of the information processed 
by the cells also increases accordingly and the rate of cor- 
relation of the cells becomes stronger. In this coupling 



scheme, 

-jr = -I- ^ kj [Pr+2j-i - Pi] 
j=o 

M 

+ Y,kr[Pr-P-i] (5) 

where, Fj — Ti{Xi,Yi,Zi) is a function and P, = 
{Xi, Yi, Zi). The extra sum indicates extra reaction chan- 
nels due to diffusively coupling of other M cells indi- 
cated by cell number index, r = 1,2,...,M, apart from 
the chain of cells. This extra term with increasing AI 
will significantly contribute to the increase of informa- 
tion transfer to the cell itself in ith position in the net- 
work topology from the coupling cells and to more other 
cells to show synchrony of more cells. 

IV. CONCLUSION 

The diffusion of inositol 1,4,5-triphosphate and cytoso- 
lic calcium ion from one cell to another through the ion- 
channels on the surfaces of the cells couple the cells giving 
rise co-ordinated behaviour of the cells. Moreover, as the 
number of diffusing molecules among the cells increases, 
the amount of information transfer among the cells is also 
increased and the cells synchronized faster. If the num- 
ber of diffusing channels is increased due to increase in 
neighbouring cells, then more cells will be coupled and 
information transfer is quicker. So we claim that the role 
of chemical coupling has significantly important role in 
the synchrony of astrocytes and long range information 
transfer among them. 
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